Excitonic exchange effects on the radiative decay time of monoexcitons and biexcitons 

in quantum dots 



Gustavo A. NarvaezQ Gabriel Bester, Alberto Franceschetti, and Alex ZungeiQ 

National Renewable Energy Laboratory, Golden, Colorado 80401 
(Dated: February 6, 2008) 

Electron-hole exchange interactions split the exciton ground state into "dark" and "bright" states. 
The dynamics of those states depends on the internal relaxation time between bright and dark states 
(spin- flip time), and on the radiative recombination time of the bright states. On the other hand, 
the calculated values of these recombination times depend not only on the treatment of correlation 
effects, but also on the accuracy of the electron and hole wavefunctions. We calculate the radiative 
decay rates for monoexcitons and biexcitons in (In,As)Ga/GaAs self-assembled and colloidal CdSe 
quantum dots from atomistic correlated wave functions. We show how the radiative decay time 
tr(X°) of the monoexciton depends on the spin-flip relaxation time between bright and dark states. 
In contrast, a biexciton has no bright-dark splitting, so the decay time of the biexciton tr(XX°) is 
insensitive to this spin-flip time. This results in ratios tr{X°) /tr(XX°) of 4 in the case of fast spin 
flip, and a ratio of 2 in the case of slow spin flip. For (In,Ga)As/GaAs, we compare our results with 
the model calculation of Wimmer et al. [M. Wimmer et al, Phys. Rev. B 73, 165305 (2006)]. When 
the same spin-flip rates are assumed, our predicted t r (X°)/t r {XX°) agrees with that of Wimmer 
et al., suggesting that our treatment of correlations is adequate to predict the ratio of monoexciton 
and biexciton radiative lifetimes. Our results agree well with experiment on self-assembled quantum 
dots when assuming slow spin flip. Conversely, for colloidal dots the agreement with experiment is 
best for fast spin flip. 

PACS numbers: 



I. INTRODUCTION: RELATION BETWEEN 
APPARENT AND MICROSCOPIC CARRIER 
DECAY 

We address here the subject of how to compare mea- 
sured exciton tji(X°) and biexciton tr(XX°) radiative 
relaxation times with calculated values. Experimentally, 
an ensemble of quantum dots is excited by an optical 
pump-pulse and the photons subsequently emitted are 
counted as a function of time. The photon emission rate 
us time is often not a simple exponential. The reason for 
this is that even in a single dot the monoexciton ground 
state is not a single state but a manifold of exchange 
and fine-structure split states with internal carrier dy- 
namics. In III-V and II- VI dots the monoexciton ground 
state originates from ej/ij, where eo and ho are, respec- 
tively, the lowest- and highest-energy confined electron 
and hole states. Due to the electron-hole exchange in- 
teraction, this state is not four fold degenerate but splits 
into four lines [Fig. nja)]iSii For the C% v symmetry of 
(In,Ga)As/GaAs self-assembled quantum dots the four 
states are the high-energy bright state B consists of a 
pair b and b' split by a few /ieV while the low-energy 
dark state D consists of a pair d and d! that is quasi- 
degenerate. B and D are split by a few hundred /xeV 
due to exchange effects. For colloidal CdSe quantum 
dots with the Cq v symmetry the internal d-d' and b-b 1 
splittings of the two pairs is small, with the b-b' recently 
measured to be about 1-2 meVi^ In turn, the b-d split- 
ting between the dark and the bright states is an order 
of magnitude larger than in self-assembled dots, ranging 
from 2-20 meV4£ In both self-assembled and colloidal 
dots, the biexciton ground state e^h^ state has no fine 



structure and corresponds to a single bright state that 
can decay to the four states of ej/ij in the monoexciton 

[Fig.mc)]. 

In this paper we show that (i) due to exchange 
and fine-structure in the monoexciton, the measured 
apparent radiative recombination time tr(X°) de- 
pends on the bright-to-dark spin-flip relaxation time 
t~bd with rate Rbd = t bd- By using an atom- 
istic pseudopotential-based approach combined with the 
configuration-interaction method^ we calculate the char- 
acteristic radiative recombination rates between B and 
the ground state (Rbo) and between D and the ground 
state (Rdo) and input them in a set of rate equa- 
tions with varying Rbd rates. We find that the pho- 
ton emission rate decays as a single-exponential with 
rate ~ Rbo for slow spin flip times; as a biexponen- 
tial for intermediate tbd\ and as a single-exponential 
with rate Rbo/ '2 for fast spin flip times, (ii) Within 
the same approach used for the monoexciton, we cal- 
culate the characteristic recombination rates Rob and 
Rod of the biexciton ground state into the bright and 
dark states of the monoexciton. We find that Rob — 
Rbo and Rod — Rdo, and that the biexciton radia- 
tive decay is a single exponential with a decay time 
t(XX°) ~ 2i?p£ regardless of Rbd- (hi) We show that 
due to the aforementioned dependence of the monoex- 
citon decay time on Rbd, the ratio tr(X°) /tr(X X°) 
has the values of 4 and 2 for the limiting cases of 
fast and slow spin flip, respectively. We thus resolve 
the apparent contradiction between the recent model 
calculations of Wimmer and co-workers^, who found 
tr{X°)/tr(XX°) ~ 2, and our previous atomistic-based 
realistic calculations of tr(X°) and tr(XX°) in which 
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we founA § t r {X°)/tr(XX°) ~ 4. We illustrate our find- 
ings with atomistic, pseudopotential-based calculations 
for a prototypical self-assembled Ino.6Gao.4As/GaAs dot 
and a CdSe colloidal dot, comparing with available data. 
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(c) Actual levels 



(d) Model levels 




FIG. 1: (Color online.) Sketch of (a) the four states that en- 
compass the monoexciton ground state ej/io and eofoo- Thin- 
and thick-line ellipses indicate, respectively, non-radiative 
thermalization as well as spin-flip channels, (b) The model 
three-level system for the monoexciton decay. Rbd (Rdb) 
is the bright-dark (dark-bright) rate while Rbo and Rdo are 
radiative decay rates for the bright and dark (model) state, 
respectively, (c) Idem (a) for the biexciton ground state eg fro 
and ej/io, and (d) model representation, with radiative decay 
rates Rob and Rod- 



II. RATE EQUATIONS FOR THE RADIATIVE 
DECAY OF THE MONOEXCITON 



Figure ^a) shows the monoexciton and biexciton en- 
ery levels that enter our calculations. We do not consider 
higher-lying states because most time-resolved photolu- 
minescence experiments are conducted at temperatures 
such that the occupation of those states is negligible. 
From Fig. ^a) we see that there are a number of dis- 
crete transitions channels denoted below by rates Rij. 
We next set up a set of channel-specific rate equations 
describing how the individual levels of Fig. ^a) "com- 
municate," from which we will deduce the global decay 
of the ground state no (i) which is measured. Using the 
characteristic radiative rates of the four excitonic states 
of ej/ij and the ground state e^/ip we establish the fol- 
lowing system of rate equations: 



dn b /dt = - {R w + R bd + R bd > + Rbo) n b + Rd'b n d > + Rdb n d + Rb'b n b > 

dn b i /dt = - (R b i b + R Vd + R vd ' + R b > ) n b > + R cVW n d > + R db > n d + R bb > n b 

dn d /dt = -(R db + R d > b > + R dd > + R d o) n d + R d > d n d > +R vd n b , + R bd n b (1) 

dn d , /dt = - (i?i4 + i?i 3 + R12 + Rw) n d > + R dd > n d + R vd > n b > + R bd > n b 

dn /dt = R d > n d > + R d0 n d + R vo n b > + R b0 n b , 



where Rij are the characteristic recombination rates from 
the level i to level j. The five-level system of rate equa- 
tions [Eq. that describe the radiative decay of X° can 
be reduced to a three-level system when (i) the thermal- 
ization rate [Fig. ^a)] within the 6-6' bright and within 
the d-d' dark states is assumed equal: 

Rbb' = Rb'b = Rdd' = Rd'd = Rth'i (2) 

(ii) spin-flip rates between the dark and the bright states 
[Fig. H[a)] are assumed to be independent on the index 



of the bright or dark state while keeping the distinction 
between bright-dark and dark-bright transition rates: 

Rbd = Rbd' = Rb'd = Rb'd' = RbD, (3) 

R d b = Rdb' = Rd'b = Rd'b' = Rdb', (4) 

and (iii) the decay rate of 6 and b' to e^hg = |0) are equal, 
and so are the decays of d and d! to |0): 

Rbo = Rvo = Rbo, Rdo = Rd'o = Rdo- (5) 
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TABLE I: Calculated values for the radiative characteristic rates Rio and Ru (i = b, b' , d, d!) for the monoexciton (X ) and 
biexciton (XX ), respectively, in an alloyed Ino.6Gao.4As/GaAs dot (base diameter b = 252 A and height ft, = 35 A) and a 
CdSe colloidal dot (diameter D = 38 A). The approximate rates Rbo, Rdo that enter the model 3-level system of rate equations 
[Eq. ©] are also shown, as well as the rates Rob and Rod for the biexciton. 

Ino.6Gao.4As/GaAs CdSe 

Rbo 0.89 ns" 1 ~ ~~ R b0 0.12 ns" 1 ~ ~~ ~ 

Ry 0.91 ns" 1 Rbo = °' 9 nS *Vo 0.12 ns- 1 Rbo = °' 12 nS 

R d0 0.65 10" 5 ns" 1 R d0 O^lO^ns- 1 „ _ 

it<j'o 0.53 10 ns fc d / 0.22 10 ns 



AA° 



#66 0.83 ns-" 1 ~ ~~ _! Ebb 0.13 ns- 1 ~ ~~ I 

D n oc -i ^os = 0.84 ns D _i i?os = 0.13 ns 

R b y 0.85 ns R bb i 0.13 ns 

i? M 0.67 10" 5 ns" 1 _ . R bd 0.09 10" 6 ns" 1 _ n 

Rbd' 0.46 10 ns it^/ 0.53 10 ns 



Assumption (i) is justified in different range of temper- 
atures determined by the magnitude of the small split- 
tings between b-b 1 and d-d': T > 2 K in self-assembled 
dots and T > 20 K in CdSe colloidal dots. Regarding 
(ii), the small fine-structure bright-dark splittings (10- 
300 /xeV) in (In,Ga)As/GaAs dots allows us to make 
the assumption Rbd = Rdb- For CdSe colloidal dots 
the bright-dark splitting is about an order of magnitude 
largeri£ and therefore Rbd = Rdb for T > 30 K. (In 
the results discussed in Sec. [V] we adopt the regime 
in which Rbd = Rdb-) Assumption (hi) is supported 
by our atomistic pseudopotential-based calculations (see 
Sec. H3). Assuming Eqs. 101-10), we simplify Eq. JTJ to 



driB 
dt 
driB 

dt 
dno 

~dt 



-{Rbo + 2Rbd) n B + 2R DB n D 
-{Rdo + 2R DB )n D + 2R B d n B 
Rbo n B + Rdo no, 



(6) 



where ns = nb+ny is the occupation of the bright states 
and riD = n^ + nd 1 the occupation of the dark states. We 
will calculate Rbo and Rdo from the electronic structure 
of the dot and vary Rbd in a wide range from very fast to 
very slow spin-flip rates to examine the different regimes 
of behavior. We will then solve Eq. © and calculate 
the phonon emission rate I(t) — Rbo nsif) + RD0 n D(t), 
which is directly comparable to time-resolved photolumi- 
nescence (PL) experiments. 



III. RATE EQUATION FOR THE RADIATIVE 
DECAY OF THE BIEXCITON 



The biexciton has a non-degenerate state without B- 
D splitting and it decays into the bright and dark states 
of the monoexciton. This decay can be modeled, simi- 
larly to the monoexciton decay, with a three-level system 
[Fig. Old)], yielding a single rate equation describing the 



population of the biexciton ground state: 
dn B {XX°) 



dt 



2(R OB + RoD)n B {XX ). (7) 



As in the monoexciton case, we will calculate Rqb and 
Rod from the electronic structure of the dot. 



IV. CALCULATION OF BRIGHT AND DARK 
RECOMBINATION RATES FROM ELECTRONIC 
STRUCTURE 

We use the empirical pseudopotential method where 
a superposition of screened atomic pseudopotentials are 
used to describe the quantum dot potential^ We take 
spin-orbit interaction into account and the method nat- 
urally includes inter-band coupling and inter- valley cou- 
pling. Following the diagonalization of the single-particle 
Hamiltonian, we use a configuration-interaction (CI) 
approach 3 to obtain correlated monoexciton and biex- 
citon wave functions |* (iy) (x)> (x = X°,XX°). The 
characteristic radiative recombination rates Rif(x) are 
calculated using Fermi's golden rule from the corre- 
lated exciton wave functions as follows. For a transition 
l* W (x)) -> \y (f) {x ~ 1)), Rifix) follows from both the 
magnitude of the dipole matrix element of the transition 
|M^(x)| and the recombination energy u>if. Namely, 



TOn C 



3 ft 2 



e— x,y,z 



(8) 

Here, e and too are the charge and mass of the elec- 
tron, respectively, and c is the velocity of light in vac- 
uum; the refractive index n of the dot material accounts 
for the material's effects on the photon emission; and 
G = G(ei n ,e ou t) accounts for the dielectric constant 
mismatch between the dot material (tin) and medium 
(e<mt) — solid barrier in self-assembled dots and liquid sol- 
vent in colloidal. 

Table [I] shows the calculated characteristic radiative 
recombination rates Rio and Ru (i = b,b',d,d') for the 
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Population of the monoexciton bright state 
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FIG. 2: (Color online.) Population n_g [Eq. logarithmic 
scale] of the model bright state of the monoexciton versus 
time (t; units of Rgg) for different of the spin-flip rate Rbd- 
Rbo = 0.9ns" 1 for an Ino.6Gao.4As/GaAs quantum dot with 
base diameter b — 252 A and height h = 35 A and Rbo = 
0.12 ns _1 for a CdSe dot with diameter D = 38 A (Table 



Rates y s and y F 




1 2 3 4 5 



FIG. 3: Fast [ 7F ; Eq. JTUJ] and slow [7s; Eq. JTTJ] com- 
ponents of the monoexciton population n_g in units of Rb- 
Inset: Amplitudes F [Eq. (O] and S [Eq. lO] versus Rbd- 
The rates Rbo and Rdo that enter 7_f and 7s are those of 
Table □ 



monoexciton and biexciton, respectively, in a prototypi- 
cal lens-shaped Ino.6Gao.4As/GaAs self-assembled quan- 
tum dot with base diameter b = 252 A and height 
h = 35 A, and a colloidal CdSe quantum dot with di- 
ameter D = 38 A. We find that the rates for bright 
states indeed satisfy Rbo — Rb>o and that dark states 
obey Rdo ~ Rd'o ~ 0. 



V. ANALYTIC SOLUTION TO THE MODEL 
RATE EQUATIONS FOR THE MONOEXCITON 

The calculated radiative rates for the 
Ino.6Gao.4As/GaAs self-assembled dot and the CdSe col- 
loidal dot (Table Pi show that it is a good approximation 
to consider Rdo = 0. In this case, together with the as- 
sumptions that Rbd = Rdb andSns(O) = nu(0) = 1/2, 
the solution of the model three-level system of rate equa- 
tions [Eq. 10] gives the following (the general framework 
for the analytic results is in Appendix |A"|> : 

n B {t) =Fexp(- 7F f) + Sexp(-7si) (9) 

with 



2 V If - is ) y ' 

S = -if ** -™ ). (13) 

In time-resolved PL experiments the measured signal I(t) 
is proportional to the number of photons per unit of time: 
dng/dt [Eq. Ip|l]. which under the assumption of -Rdo = 
results in I(t) = R B oti B (t). Note that by being propor- 
tional to the occupation of the bright state the signal I(t) 
carries information on both radiative and non-radiative 
(spin flip) processes. Figure |21 shows the logarithm of 
nsit) as a function of time for different spin-flip rates 
Rbd, and Fig. |31shows the slow (75) and fast (7^) com- 
ponents of the decay of n B (t) versus Rbd- We find that 
in the limiting cases of (1) extremely slow {Rbd <C Rbo) 
and (2) extremely fast (Rbd ^> Rbo) spin flip the decay 
of tib is primarily determined by a single exponential. In 
case (1) we find 7f — Rbo + ^Rbd and 7s ~ 2Rbd, while 
F 1/2 and S ~ [Fig. [3J inset]. Here, jf resembles 
the expected result for the decay rate of the PL in the 
presence of nonradiative recombination centers; the dark 
state of the monoexciton in this caseii2iii The population 
of the bright state is 



If = \{Rbo + ^Rbd) + \^R 2 bo + ( ar bd) 2 , (10) n B (t)~-exp[-(R B0 +2R BD )t]. (14) 

is = -^(Rbo + 4:Rdb) - 2\/ r bo + {^ r bd) 2 (11) i n this regime, I(t) decays approximately with the char- 
acteristic lifetime of the bright states. In case (2), we find 
and 7f — Rbo/2 + ARbd and 75 ~ Rbo/2, and S ~ 1/2 and 
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F ~ [Fig. |2| inset]; therefore, 



VII. THE RATIO tr{X°)/tr{XX°) 



u\ 1 ( R- m l 
n B [t) ^ -exp I — i 



(15) 



thus decays with an approximate characteristic time 

of T R {X°) 

lifetime of the bright state 



2R B ^; twice as large as the characteristic 



VI. COMPARISON WITH EXISTING 
EXPERIMENTAL DATA 

In the experimental literature, data on the bright-dark 
transition time are scarce. (In,Ga)As/GaAs dots. Dal- 
garno and co-workers*^ have studied recently the effect 
of the dark state in the decay of the monoexciton in a 
gated structure and have estimated Rbd > 20 ns at a 
temperature of 5 K. In addition, they have found that 
this spin- flip rate varies strongly with the applied bias. 
From time-resolved PL experiments, Favero et ali& have 
extracted significantly disparate values for two different 
dots: R B d ~ 440 and 30 ns. 

InP/(In, Ga)P dots. In a two-photon absorption exper- 
iment, Snoke and co-workers^ populated the dark states 
of the monoexciton and measured the subsequent lumi- 
nescence as a function of time. From the data below 
70 K, the authors found that the spin-flip, bright-dark 
transition time satisfies Rbd — 200 ps. 

CdSe dots. By performing fluorescence transient ex- 
periments at room temperature, Wang et alm& have con- 
cluded that the bright-dark, spin-flip relaxation times 
R 



BD 



0.2-0.4 ps. Thus, the spin-flip process is or- 
ders of magnitude faster than the characteristic radiative 
recombination time. These experiments reveal order-of- 
magnitude variations in the bright-dark spin-flip times, 
suggesting that the value of Rbd appropriate to inter- 
prete time-resolved photoluminescence experiments in 
quantum dots is controversial and further research is 
needed to understand the spin-flip mechanism. 

In the slow spin flip regime, consistent with the find- 
ings of Dalgarno et al. 9 and Favero et al. 13 , our calculated 
value t r {X°) ~ R B q = 1.1 ns for the Ino.6Gao.4As/GaAs 
dot is in excellent agreement with the data of Bardot et 
a/Ji who extracted 1.55 ns from time-resoved photolu- 
minescence, and the value of 1 ns found by Buckle et 
aliU and Stevenson et ali& In the fast spin flip regime, 
for our prototypical CdSe dot (see Table P) , we obtain 
t r {X°) — 17 ns, which is in excellent agreement with 
the values of 17 ns and 19 ns extracted, respectively, by 
Brokmann et al. 1 ?- and Labeau et alm& from time-resolved 
photoluminescence in ZnS-passivated CdSe dots. In the 
regime of intermediate spin flip rates, measurements of 
the biexponential decay of I{t) — as those performed by 
Dalgarno et al£ — could be used to deduce the spin-flip 
rate. 



To compare the biexciton decay rate with to the mo- 
noexciton decay rate in the limiting cases discussed above 
(Sec. EJ, we first note that atomistic pseudopotential- 
based calculations show (Table P) that the characteristic 
radiative rates for the biexciton ground state satisfy 



Rqb(XX°) ~ Rbq{X ), 
Rod{XX°) ~ Rd (X°). 



(16) 



Second, in contrast to X°, we note that the solution 
of the rate equation for XX [Eq. J7J] results in 
riB(XX ) ~ exp(— 7t); a single exponential that decays 
with rate 



7 = 2(R 0B + Rod) 2(R B0 + R D0 ) 



(17) 



regardless of the value of the spin- flip rate Rbd- Similarly 
to X°, the time-resolved PL signal is proportional to the 
population of the bright state of the biexciton. For slow 
spin flip [case (1), Sec. |V] we find a decay-rate ratio 
between X° and XX of 



2Rob 



2 (slow spin flip), (18) 



Rbq + 2Rbd 
and for fast spin flip [case (2), Sec. |V] we find 



7/7S 



2Rqb 
Rbq 1 2 



4 (fast spin flip). 



(19) 



We emphasize that depending on the magnitude of the 
spin-flip time the ratio t r {X°) / t r {X X°) can change by a 
factor of two, therefore the assumed spin-flip time is cru- 
cial when comparing results for tr(X )/tr(XX ). Re- 
cently, Wimmer et al£ have used a quantum Monte Carlo 
(QMC) approach with model single-band effective-mass 
electron and hole states to calculate 



t r (X°)/t r (XX°) ~ 2 {slow spin flip). (20) 

Those authors speculated that the disagreement with the 
pseudopotential and CI calculations of Ref. H, which 
adopt the fast spin-flip regime and predict 



t r (X°)/t r {XX )~4 {fast spin flip), (21) 

originates from an inaccurate treatment of correlations in 
CI. However, as is obvious from Eqs. 118|) and 119|) . the 
discrepancy can be directly attributed to the different as- 
sumptions for the spin flip rates. When the same spin-flip 
rates are assumed, our results for t r {X°)/t r {XX°) are 
in agreement with the QMC results of Wimmer et ali, 
suggesting that our treatment of correlations is adequate 
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to predict the ratio of monoexciton and biexciton radia- 
tive lifetimes. On the other hand, the calculated values of 
the radiative recombination times depend not only on the 
treatment of correlation effects, but also on the accuracy 
of the electron and hole wavefunctions. Our atomistic 
results are in good agreement with experiment while the 
results of Wimmer et al£ based on the single-band effec- 
tive mass approximation differ from experimental data 
by a factor of two. 

Finally, note that in our calculation of 
tr(X )/tr(XX°) we assume that the the change 
in occupation of the monoexciton bright and dark states 
[Eq. |JTJ] is not affected by the decay of the biexciton 
state. 



VIII. SUMMARY 

We calculated the characteristic radiative recombina- 
tion rates for the ground state of the monoexciton and 
biexciton in self-assembled (In,Ga)As/GaAs and colloidal 
CdSe quantum dots using atomistic wave functions. For 
the monoexciton we used these rates in a model three- 
level system of rate equations where we varied the spin- 
flip rate Rbd- The latter affects significantly the radia- 
tive decay time: Fast spin flip leads to an exciton radia- 
tive recombination rate twice as fast as the rate obtained 
from slow spin flip. The radiative decay times tr(X°) 
calculated in the limit of slow spin flip are in excellent 
agreement with available data for self-assembled dots, 
while for colloidal dots the agreement is best for fast spin 



J 



flip. The biexciton radiative decay is a single exponential 
with a relaxation time that is independent of the spin- 
flip rate. But the ratio between the radiative decay time 
of the biexciton tr{XX°) and monoexciton does depend 
on Rbd and results, respectively, in trX° /tr{X X°) ~ 4 
and 2 for fast and slow spin flip. This result resolved 
the apparent contradiction between the calculation of 
Wimmer et al£, who predicted trX° / tr(X X°) ~ 2 and 
our previous atomistic calculation^ in which we found 
t r (X°)/t r (XX°) ~ 4. 
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APPENDIX A: SOLUTION OF THE 
THREE-LEVEL MODEL RATE EQUATIONS FOR 
THE RADIATIVE DECAY OF X° 

We use the simplified model of Eq. JSJ as shown in 
Fig. H£b). To find nsif) and no if), we propose 



n B it) = F exp(- 7F t) + S exp(- 7s t) , (Al) 
with rifl(O) = F+S and 7 _p 75, and solve first for nof) 
with initial condition n£>(0), where ns(0) + tid{0) = 1- 



Then, to solve for S, F, 75, and 7 _p, we substitute the solution of nof) in the rate equation for dns/dt [Eq. @], 
obtaining the following conditions. 



7^ — {Rbo + 2Rbd + Rdo + 2Rob) 7f + RboRdo + 2RdbRbo + 2RbdRdo — 0, 
75 ~ (-Rbo + 2Rbd + Rdo + 2Rdb) 7s + RboRdo + 2RbbRbo + ZRbdRdo = 0, 

2R DB MO) + iR D0 + 2R DB ) n B iO) - {RboRdo + ZRdbRbo + 2R bd Rdo) ( — + — )= 0- 

\1F IS J 



(A2) 
(A3) 

(A4) 



Thus, we find 



1 



"if 



7s 



{Rbo + Rdo + 2{Rbd + Rdb)) 



+ \V{Rbo - Rdo) 2 + MRbo - Rdo)(Rbd - Rdb) + 4{R BD + Rdb) 2 
— {Rbo + Rdo + 2{Rbd + Rdb)) ~ 

-IViRso - Rdo) 2 + MRbo - Rdo){Rbd - Rdb) + ±{Rbd + Rdb) 2 



(A5) 



(A6) 



and 
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Rbo + 2R B d - 7s , n s 2Rdb , . 

F = n B {0) n D (0) (A7) 

7f - 75 If — 7s 

Rbo + 2Rbd - If , n s 2i?£, B 

6 = ns(0)H n D (0) (A8) 

7f — 7s 7f - 7s 
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